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ABSTRACT 

We investigate the three-neutrino mixing scheme for solving the atmospheric and 
LSND anomahes. We find the region in the parameter space that provides a good fit 
to the LSND and the SK atmospheric data, taking into account the CHOOZ constraint. 
We demonstrate that the goodness of this fit is comparable to that of the conventional 
fit to the solar and atmospheric data. Large values of the LSND angle are favoured and 
sin^(2^LSND) can be as high as O.L This can have important effects on the atmospheric 
electron neutrino ratios as well as on down-going multi-GeV muon neutrino ratios. We 
examine the possibility of distinguishing this scheme from the conventional one at the 
long baseline experiments. We find that the number of electron neutrino events observed 
at the CERN to Gran Sasso experiment may lead us to identify the scheme, and hence 
the mass pattern of neutrinos. 
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1 Introduction 



The present data from the experiments on atmospheric, solar and accelerator (LSND) 
neutrinos indicate neutrino flavour oscillations. The data from each of these sets of ex- 
periments individually can be explained by a single dominant mass square difference Am? 
and a mixing angle 6 between two active neutrinos. The atmospheric neutrino data from 
SuperKamiokande (SK) [|I| indicate z/^ ^ Ur as the dominant mode, with Am^^^^ = (1 
- 8) X 10-3 eV^sin2 2^ = 0.8 - 1.0. The three MSW solutions (LMA, SMA and LOW) 
as well as the vacuum oscillations can provide reasonable fits to the solar neutrino data 
[|, |, |, I and all these solutions have Am| < 2 x 10"^ eV^. The results of the LSND 
experiment [^, are neither confirmed nor fully excluded by the KARMEN2 data |^ , and 
the combined fit allows a region of Am^^^^ = (0.1 - 1) eV^,sin^2^^e ~ lO'^ — 10"^. 



There are also constraints on the mixing of Ue from the CHOOZ experiment [|10|: we have 
sin2 2^e < 0.1 for Am^ > lO'^ eV^. 

In the context of only three known neutrino species, the Am^s corresponding to the 
solutions of the three neutrino anomalies above (atmospheric, solar and LSND) cannot be 
reconciled. The three-neutrino schemes make a very poor fit to all the data. Two ways 
are possible out of this predicament: (i) turn a blind eye to one of the experiments and fit 
for the other two in three-neutrino schemes; (ii) solve two of the neutrino anomalies with 
three-neutrino oscillations and solve the third one by using exotic models such as sterile 
neutrinos, FCNC, neutrino decay, extra dimensions, etc. 

Since the LSND result is yet to be confirmed, it is customary to ignore it and accom- 
modate the solar neutrino deficit and the atmospheric anomaly with the mixing between 
three active neutrinos. Implicit in this is the assumption that the LSND results will be 
proved false by future experiments, which can be justified by the fact that KARMEN2 
[§] and Bugey already rule out most of the allowed region of LSND. However, this 
is just a convenient assumption, and the possibility of the LSND results being confirmed 
by future experiments such as BooNe cannot be ignored. Also the new analysis of 
the final data by the LSND collaboration |^ is consistent with the old results [^], and 
therefore strengthens the anomaly evidence. 

Our approach will be to study the neutrino anomalies with three-(active)-neutrino 
oscillation. This allows us to solve only two of the three anomalies. As the atmospheric 
data are showing strong evidence for neutrino oscillation, thanks to the large range of 
L/E probed, we will take this to be one of the anomalies solved by oscillation. There is 
no compelling evidence that the electron neutrinos participate in the oscillations of atmo- 
spheric neutrinos. This implies that the P^^^y^ must be small, meaning either that the 
mixing angle is small (LSND case) or that the Am^ is too small to affect the atmospheric 
neutrinos (solar case). The large angle solutions to LSND is in any case ruled out by the 
results of Bugey. Moreover, just from the point of view of goodness of fit (quantified by 
a function), the best fits to (I) atmospheric and solar data, and (II) atmospheric and 
LSND data are equally good (as we shall show in this paper). Scheme II gives a differ- 
ent mass spectrum from the conventional scheme I, and an eye should be kept on future 
experiments in order to resolve this discrete ambiguity. That we leave out one anomaly 
does not mean that we do not believe in it, but rather that it has to be solved in some 
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other way. 

In the case of solar neutrinos there are various viable "exotic" solutions. The most 
popular among them is the sterile neutrino SMA MSW solution. In the favoured "2+2" 



four- neutrino scheme ||T3[, the sterile neutrino participates mainly in the solar neutrino 



anomaly and the three active neutrinos solve the atmospheric and LSND anomalies among 



themselves |15]. This is indeed a possible extension of scheme II, and the fit we perform 



will be a guide for this option. Many other exotic solutions exist for solving the solar 



neutrino anomaly |14]. Flavour changing neutral current interactions (FCNC) [0 can 



explain the data well |T^ - this solution uses the matter of the Sun to make a FCNC 
transition, such as v^. + quark quark + v^. The different neutrino production points 
gives the necessary energy dependence of the solar neutrino flux. Furthermore a violation 
of the equivalence principle has been suggested and cannot be excluded. Solutions have 
also been suggested in the scenarios with large extra dimensions, where the solar neutrino 
anomaly is accounted for by the mixing with a Kaluza-Klein tower of sterile neutrinos 
and the three active neutrinos account for the atmospheric and LSND anomalies 



r^]. A resonant spin- flip conversion pO] of the left-handed neutrino to unobservable 



right-handed states, due to the solar magnetic field, gives a good fit to the data [21 



This solution cannot be reconciled with LSND and atmospheric anomalies, as the mass 
squared difference needed is too small (~ 10~^eV^). In the three-neutrino schemes able 
to account for the atmospheric and LSND result, one of the possible mass patterns has 
the electron neutrino mainly in the heavy state, and therefore neutrino decay could also 
be thought of as a way out for explaining the solar anomaly. Although this solution does 
not agree well with the data. 

Some exotic solutions for the LSND anomaly have also been suggested in the literature 
for instance, having new flavour- violating decay modes of muons. This possibility 
cannot yet be ruled out in a model-independent way. The introduction of a sterile neutrino 
to solve the LSND anomaly gives rise to the "3+1" neutrino mass pattern PB[, which is 



very close to being excluded It might be argued that the exotic solutions are more 
likely to solve the solar anomaly than the LSND one, since the matter and the magnetic 
field, which provide more degrees of freedom for the exotic solutions in the solar case, are 
absent in the LSND experiment. 

There is a vast literature on three- neutrino oscillation fits to atmospheric data [ p5|] . 
Also the combined fit to solar and atmospheric data is well-explored and is known to 
describe the data well [^, |2^. Furthermore the two- neutrino z/^ — > oscillation gives a 
good explanation of the atmospheric data. This solution would be obtained in the limit of 
a very small LSND angle, and therefore we already know that the LSND and atmospheric 
data can be reconciled. The important question to answer is if the three-neutrino scheme 
can fit the data even better and if there are any possible effects that could be observed 
in the near future. 

The three-active-neutrino solution to LSND and atmospheric problems has also been 
studied in ||2^, where however only the up/down asymmetries from the atmospheric data 
were used for the fits. Here we will consider the full set of 40 data points from the sub-GeV 
and multi-GeV neutrinos observed by the SK collaboration. Furthermore we take into 
account the data from CHOOZ [|l^ and the matter effects inside the Earth. We define 
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a function and perform fits to schemes I and II above, and find that both fits are 
equally good. We explore further the fit to scheme II. The mass pattern corresponding 
to this fit may be tested at the long baseline experiments, e.g. K2K ||2^, MINOS |^ 
or CNGS [^. We perform Monte Carlo simulations to check if this mass pattern can 
be distinguished from the conventional one at these experiments. If BooNE confirms the 
LSND results, this scheme II, coupled with the appropriate exotic solution for the solar 
neutrino anomaly, will provide the solution for the mass spectrum of neutrinos. 

The paper is organized as follows. In Sec. ^ we perform fits to the data for schemes 
I and II. In Sec. ^, we compare the signals at the long baseline experiments with the 
two mass patterns and look for ways of distinguishing between them. Our findings are 
summarized in Sec. H. 



2 fits to the data 

The aim of this section is to find the regions in the three- neutrino mixing parameter space 
that describe (I) the atmospheric and solar data and (II) the atmospheric and LSND data. 
We shall define a function in order to quantify the goodness of these fits. It is clearly 
not fair to compare the x^s from these two fits, since they correspond to different data 
sets. We perform the fit to scheme I to demonstrate the soundness of our x^ function by 
showing that this procedure, rough though it is, reproduces the conventional fits, which 
are much more accurate Since our aim here is just to get a broad estimate of the 

allowed parameter space, this rough fit should suffice for our purpose. Then we follow the 
same procedure to perform a fit to scheme II and find the allowed region in the parameter 
space that describes the data reasonably well. 



The neutrino mixing matrix Umns is parametrized as |32 



U. 



MNS 



' C12C13 S12C13 Si3 \ 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 (1) 
\ S12S23 — C12C23S13 — C12S23 — S12C23S13 C23C13 j 



where we have neglected any CP violation. The mass spectrum is parametrized by five 
parameters, S12, S23, S13, l^m^^ and lS.m\^. (Here and in the following Sj^, Cjj are short- 
hands for sin(%), cos(%) respectively, and we use the convention Am^^ = — m^.) We 
define v\ as the lightest state and z/3 as the heaviest state. 

For both schemes there are two solutions; the normal hierarchy (Amgi -C Am^g — 
Arrig^) and the inverted hierarchy (Amgj^ ~ Am^^^ ^ I\vn^^. Neutrino oscillations in 
vacuum cannot distinguish between the two hierarchies, but with matter effects, the data 
could in principle distinguish between them. In scheme I the characteristic features of 
the normal hierarchy is the small value of S13, and the inverted hierarchy corresponds 
to values of S12 close to 1. It has recently been shown that the differences between the 



two hierarchies within this scheme are very small ||33|. For scheme II normal (inverted) 
hierarchy corresponds to S13 ~ 1 (si3, S12 — 0). The matter effects in this scheme are 
negligible, since both the Am^^^^ and AmLg^D too large for the Earth's density to 
play any significant part. Therefore, the differences between the predictions from the 
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two hierarchies are extremely small. For both schemes we shall only consider the normal 
hierarchy while performing the fit. It should be remembered that for each point in the 
parameter space with normal hierarchy, there exists a corresponding point with inverted 
hierarchy that gives almost the same value of x^- 

For the atmospheric data, we only consider the SK multi-GeV and sub-GeV data 



3J]. The mean neutrino energy for the partially contained events and the upward going 
muons at SK is higher. In this range the effects of a solar mass squared difference are 
therefore negligible. The small contributions arising from an averaging of a LSND mass 
sqaured difference can be compensated by a shift towards smaller Am^tm (see section pT^ ) . 
Therefore it is not expected that the data will affect the comparison of the two schemes. 
Nevertheless it might results in small changes in the allowed regions. The further inclusion 



of the data from other atmospheric neutrino experiments |^ is not expected to affect 
the fit much, since the number of fully contained events at SK is overwhelmingly large as 
compared to these. 

The experimental data are represented by the ratios Ra^i between the experimental 
values for the fluxes and the theoretical Monte Carlo prediction in the case of no oscillation 
for muon and electron neutrinos in the 10 different zenith angle bins. The ratios Ra,i can 
be written as 

RaI = K7/K^ , a = /.,e , z = 1...10. (2) 
The number of events for the sub-GeV neutrinos is calculated as; 

E poo rl rEi.rmi^ nl I p2iv 

> / dEy I dcos{9u) / dEi / c/cos(6';^)— / dcj) 

^ ,J^T'(,^,J^f'''(, , Pc^p(EM<Ei) , (3) 
dEydcos[Oy) dEidcos{Oi^) 

where Ei, is the neutrino energy, Ei is the lepton energy, 6iu is the angle between the 
neutrino and the scattered lepton, is the zenith angle of neutrino. The number of 
target nucleons is denoted by ut- <p is the azimuthal angle of the incoming neutrino, and 
is used to calculate the zenith angle 6i of the charged lepton: 

cos{9i) = cos{6^) cos{6^i) + sin(^^) cos(0) sin(6';) . (4) 

For the differential sub-GeV cross section for uiN IX, we use the quasi-elastic approx- 
imation [^, with the sub-GeV fluxes $ taken from [^. The efficiency function e{Ei) is 



not published by the SK collaboration. But as we divide by the SK Monte Carlos, this 
effect is averaged out if the efficiency function is flat within the different samples. 

In the calculations of the multi-GeV ratios we find that the number of events can be 
well described by 

Ni=nT dE, / dcos{e,)-—- a^{E,)P^p{E,, 9,) . (5) 

Je,„,,„ J -I dE, A cosiO,,) 



Here we assume that the energy of the charged lepton is the same as that of the incoming 
neutrino. In order to account for the scattering angle, we further smear the spectrum 
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with a Gaussian function, where we use the width a ^ 25° for the electron neutrinos 



and cr ^ 15° for the muons neutrinos |38|. The total cross section for neutrinos is given 
in [^], and for the antineutrinos we use the approximate formula (777^ = 0.34 x 10~^^Ei, 
cm^/GeV. The multi-GeV fluxes are taken from [^]. A comparison between the ratios 
obtained from and the ones obtained by using the differential cross section including 
the quasi-elastic, one pion and deep inelastic channels yields hardly any difference. 

The probabilities are calculated taking into account all the mass-squared differences: 



Pal3 = - 4 ^ U a,iU pjU ^jU pi SIV? ) , (6) 



where all the quantities are calculated in the presence of matter wherever appropriate. 
Even in the conventional case, the sub-GeV ratios in particular are influenced also by the 



small mass squared difference The oscillation length L depends on the zenith angle 
of the neutrino: 



L = ^J{Re + hf - Rl sin\e,) - Re cos(e,) , (7) 

where Re is the radius of the Earth and h is the production height of the neutrinos in 
the atmosphere. We take /i ~ 10 km. 

Matter effects are important for the atmospheric neutrinos. We simulate the matter 
effects by using a two-shell model of the matter densities in Earth. The density in the 
mantle (core) is taken to be roughly 3.35 (8.44) g/cm^, and the core radius is taken to be 
2887 km. From this we calculate the average density as a function of the neutrino zenith 
angle. This allows us to use the three-neutrino mixing matrix in matter, as calculated in 

We define the atmospheric ^is 



10 (V^^P _ Dth \2 
Xatm / J / J / J 2 



M,S Q=e,li 1=1 ^ai 

where the statistical errors and M, S stand for the multi-GeV and sub-GeV data 

respectively. The total number of data points is 40. 

From the CHOOZ experiment [|10| we use the 15 data points with the statistical errors. 
The CHOOZ baseline is roughly 1000 m, and the neutrino energy in the range from 3 to 9 
MeV. Therefore the sin^(Am^L/4ii^) terms do oscillate within the energy region and hence 
all the data points need to be taken into account. In each bin we average the probability 
over energy. When including the LSND experiment we use one datum p. The distance 
travelled by the antineutrinos (z7^) is set to 30 m and we use the mean energy of 42 MeV. 



When including data from the Bugey experiment |]1T| we take three data points; P^eFe for 
the three different baselines of 15 m, 40 m, and 95 m. The probability is averaged over 
the positron energy range 1-6 MeV. We take both the statistical and systematic errors, 
as the systematic errors are large compared to the statistical ones. 

The solution to the solar neutrino problem is preferred to be within the LMA region 
and we will therefore use parameters only within this region in order to make a fit. The 

for the solar LMA solution is taken from Ref. ||26||. This fit includes the total rates 



in the chlorine (Homestake) experiment [^, the gallium (SAGE, GALLEX) experiments 
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1^, Kamiokande |^3[ and SuperKamiokande [Q]. Furthermore the spectra for both day 
and night are included. In total there are 42 data points fitting with three parameters: 
Amg, 9q and the overall flux normalization for the SK spectrum. 
The total function for scheme I is then 

Xl — Xatm + XcHOOZ + XhMA ; wJ 

and for scheme II, we have 

Xjl = Xatm + XCHOOZ + XlSND • (10) 

Before interpreting the results, a word of caution is in order. In our analysis (as 
in the standard procedures), we deal not with the number of observed events but with 
the ratio of the number of observed events to the number of events expected from the 
Monte Carlo. These ratios are convenient because they directly give a measure of the 
survival/oscillation probability. The Monte Carlo predictions, in the case of atmospheric 
neutrinos fluxes, contain large uncertainties, especially in the absolute values of fluxes. 
This is the reason why some data are presented in the form of a ratio of two measured 
quantities and is compared with theoretical predictions of the same ratio. We have not 
included theoretical uncertainties and correlations in the atmospheric neutrino fluxes. 
However, a previous analysis has found that they do not affect the fit significantly. 



It has been pointed out recently that there is a significant discrepancy between the 
commonly used primary cosmic ray fluxes and the measured ones |^^. The variation of 



the primary cosmic ray flux is directly related to the absolute value of the atmospheric 
neutrino fluxes and can also affect (to a less significant extent) its energy and zenith angle 
dependence. This may affect the allowed region of the atmospheric neutrino anomaly. It 
is important to bear in mind that when applying the calculated atmospheric neutrino 
fluxes to a neutrino oscillation study, the absolute values of the fluxes as well as the tiny 
variations with energy between them could become important. 

2.1 Fit to solar, atmospheric and CHOOZ data (scheme I) 

In this section, we perform a fit to the conventional three-neutrino scheme for solving the 
solar and atmospheric anomalies. When fitting within this scheme we will take Amg^ ~ 
10~^ - 10^^ eV^ and lS.m\^ suitable to solve the solar neutrino problem. The Xatm 
XcHooz S'^s to a very large extent independent of S12. We will therefore keep the solar 
angle constant at — 0-3; which is close to the best-fit point of the LMA solution 
The small ^m^^ is taken in the range 2 x 10~^ - 2 x 10~'^ eV^. In particular the 



electron neutrino ratios can be affected by the solar mass difference |^T|, and we therefore 
calculate Xatm + Xchooz ^ ^w values of Ami,. We thus have four degrees of freedom, 
S23,si3, Am^i, Am^P 

We first find the best point while doing a combined fit to SK and CHOOZ, with 
= xltm + XcHooz- There are 55 data points: 40 from SK atmospheric data and 15 
from CHOOZ. The minimum is Xmm = 45.3 at^ 

= 0.41 , = 0.00 , Am^i = 2.0 x 10^^ , Am^^ = 4.1 x 10"^ . (11) 



^ Henceforth, we implicitly assume the units eV^ for Am^, unless specified explicitly. 
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Note that the minimum occurs for the largest value in the allowed region for AitIq. With 
Am^i constrained to be 2.0 x 10"^ (8.0 x 10"^), the value of Xmm is Xmm = 49.1 (46.9) at 

= 0.58 (0.40), s?3 = 0.01 (0.00), Am^^ = 0.0040 (0.0045). We see that the dependence 
on the solar mass squared difference is weak. The rises slightly when Am^ becomes 
smaller. For AttIq < 10~^ this dependence is lost and the function is almost constant. 

Figure |l|(a) shows the allowed region in Am|^ = Aml^^ and iU'^^U^^, ~ sin^(26'atm)- 
The confidence intervals are calculated for four parameters (as we keep Si2 fixed), so that 
the 90% (99%) confidence interval corresponds to - Xmm < 7.8 (13.3) |1. At 99% CL 
we find 1.75 x lO^^gy^ < Am^ < 7.5 x lO'^ eV^ and sin2(2^atm) > 0.83. These bounds 



are in good agreement with previous analyses [25 



In figure |l|(b) we show the 99% CL and the 90% CL area in sj^ and sj^. At 99% CL the 
upper bound on sf^, is 0.065, corresponding to a value of sin^(26'cHooz) = 4:11^^(1 — U^^) 
smaller than 0.25. 

The full three-neutrino fit is not symmetric around = 0.5, as is also seen in Fig. H(a) 
and (b). The small asymmetries arise from non-zero values of either S13 or AttIq. This 



may be interpreted in terms of the electron excess (see the discussion in Sec. ^ 



We now consider the combined fit including the effect of the solar neutrino data in 
the LMA region. The common parameter, AitLq, must be taken to be the same, and we 
calculate xj for a few values of this parameter. The minimum value occurs at AitLq = 8 x 

10~^ with Xlma ~ 37.3 Hence when fitting to all three experiments (SK atmospheric, 

2 



CHOOZ, solar data), we get X/(min) ~ ^4 with 6 parameters and 97 data points, so that 
X?(mm)/dof ^ 0.93. 

Recently a fit was performed to the solar data allowing for a free B and hep flux in 
the Sun The obtained LMA x^ was 29.0 for 39 data points with 4 fit parameters 



AitLq, 6q and the fluxes of and hep neutrinos. Combining this with our fit for SK and 
CHOOZ, we would get X/(min) ~ 78 , so that X/(min)/'^of ~ 0.90. Allowing for an overall 
normalization of the atmospheric fluxes, as done by the SK collaboration, could in general 
also improve the goodness of fit for atmospheric data. Let us also note that in Ref. |^| 
a combined fit to atmospheric, solar and CHOOZ data gave a best fit oiQ x^/dof = 0.97. 

The LMA solution currently gives the best fit to the solar data, and the fit to the 
atmospheric data is not affected much by very small solar mass differences. Therefore we 
expect that the value of x^ can only be higher in other regions of parameter space than 
the one computed above. 



2.2 Fit to atmospheric, LSND and CHOOZ data (scheme II) 

In this section we give the results of fitting to the SK atmospheric, CHOOZ and LSND 
data and hence will take Amli = Aml^^ and Am|^ ~ Amfj = ^^Ism- The SK and 
CHOOZ data are largely independent of the large Am^, which we can therefore freely 
use to fit the LSND experiment. For every point in parameter space we simply choose to 
calculate the value of Amg^^ that fits the LSND data best. The remaining four parameters 
(si2, S23, S13, Am2i) are varied to obtain confidence levels. 

^The difference from our fit is that in j26j, the Amg effects on xitm ^''^^ neglected, whereas we neglect 
si3 effects on Xlma- 
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Let us first perform a combined fit to the SK and CHOOZ data with = xltm + 
XcHooz' we did for scheme I in Sec. |2?1| (again we have 55 data points). The confidence 



levels are calculated for four fitting parameters. Fig.^(a) shows the allowed region in 
Aml^^ and 4:U^iU^2 - sm^{29^tm)- At 99% confidence level the bounds are 4:U^iU'^2 > 0-8 
and 1 eV^ < Aml^^ < 7 eV^. The allowed region is very similar to the one found in the 
conventional case. It shows that the two-neutrino atmospheric neutrino fit is only slightly 
affected by the three-neutrino extension in both schemes. Although we note that the 
region in Aml^^ is lowered in scheme II. This is due to the small constant contributions 
that arise from the averaging of the sin^ {Amlg^-^L / 4:E) terms. 

Let us now include the LSND data in the fit. The best-fit point with ^//(min) ~ ^^-^ 
for the combined fit is obtained at 

s^2 = 0.69 , = 0.98 , = 0.98 , Amj^ = 0.0031 , Amj^ = 0.22 . (12) 

This should be compared with the two-neutrino z/^ — > Ut minimum value of of 49.5 
(at Aml^^ = 0.0038). We remind the reader that the LSND datum is chosen to be fitted 
best. Hence the effects of the three-neutrino scheme are not large, but still significant. 
The main reason for the lowering of the is an excess in both sub-GeV and multi-GeV 
electron neutrino ratios. The goodness of fit to the data, having Xjj/dof ~ 0.88 is thus as 



good as the one obtained to the solar, atmospheric and CHOOZ data in Sec. |27l . 
The LSND probability is 

P,,^., ^ AU!,Ui, sin^ l^^ia^^ _ AU,,U.2U,,U,2 sin^ (^^) , (13) 

where the second term is small because of the small value of sin^ {Aml^^L / AE) , so that 
we may define the LSND angle as sin^(26'LSND) = ^U^sU'^s- In terms of physical angles, 
the best-fit point corresponds to sin^(26'atm) = 0.99 and sin^(26'LSND) ~ 0.08. 



The best-fit point ([13|) has an LSND angle in the region excluded by Bugey ||rT[. We 
investigate the impact of the Bugey data on the value of the LSND angle when combining 
these in the fit. Including the Bugey data gives a best-fit point with ^//(min) — 47.0 at 

= 0.7 , = 0.98 , = 0.985 , Am^^ = 0.0035 , Aml^ = 0.26 . (14) 

or in terms of physical angles sin^(26'LSND) — 0.06, sin^(26'atm) — 0.98. With 59 data 
points fitted with 5 parameters, we have x^/dof ~ 0.87. The upper bound at 90% CL is 
sin^(26'LSND) < 0.11. Therefore the inclusion of the Bugey data still favours large LSND 
angles. 



Also the best-fit point (O) is slightly different from the one obtained in |28|; in partic- 



ular the LSND angle is larger. It should be said that in Ref. the LSND data were not 
included in the fit, but only a mass squared difference of order eV^. We therefore believe 
that the discrepancy is due to the restriction Am^^^y^ > 0.5 in |]2S|, which is equivalent 



to restricting the value of the LSND angle to sin^(26'LSND) < 0.02, and also to the use of 
up/down asymmetries as representing atmospheric data. 

Since Ues = sis ~ 1 in scheme II, the solutions are still close to an effective two- 
neutrino mixing between i/^ and Ur, with dominating the third mass eigenstate. For 



8 



the normal hierarchy (discussed here), Vf, is almost equal to the heaviest state, whereas 
for the inverted hierarchy, Vf. is almost equal to the lightest state. 

Another interesting point about scheme II is the connection between the "LSND angle" 
and the "CHOOZ angle" . The LSND angle has already been defined above. In order to 
define the "CHOOZ angle", let us have a look at the CHOOZ survival probability for 
Vf>. Since the matter effects are small, for the analytical approximations we may put the 
parameters equal to their vacuum values. The z/g survival probability can be approximated 
by 

P..^.. = 1 - mi - U!,) sin^ ^ ^^'^SNoL ^^ _ (^^) 

^ l-2UUl-U!,) . (15) 

The dominant term in the survival probability Pge originates from the Attt-lsnd term, 
since the coefficient in front of the Aml^^ is very small (of order cf^). Hence the CHOOZ 
angle is approximated by sin^(26'cHooz) ~ ^U^si^ ~ U^s)- 

Figure ^(b) shows the allowed region in the parameter space of sin^(26'cHooz) and 
sin^(2^LSND)- The lower (upper) bound on the LSND angle then turns into a lower 
(upper) bound on the CHOOZ angle. When considering the combined fit to the atmo- 
spheric, CHOOZ and LSND data, the upper bound to 90% CL with 5 parameters is 
sin^(26'cHooz) < 0.20. The corresponding lower bound is sin^(26'cHooz) > 1-2 x 10~^. 
Although the bound is very weak, it does require in principle a full three-neutrino mixing 
scheme, and will be important in future experiments trying to extract CP-violation. 



Let us briefiy note that the scheme suggested in |^^, which attempts to solve all 
three neutrino anomalies with just three fiavours, is found to be strongly disfavoured, in 
agreement with previous results |H8[. Within this scheme, is always above 120. Recently 



the authors of |49| also suggested two different regions, which they found to be good 



candidates for solving all three anomalies. We find these regions strongly disfavoured; in 



particular the region suggested in Table 2 of confiicts with the CHOOZ data and can 



be considered ruled out. The of the SLMA solution defined in is also too large 
to be consistent with the data: at = 0.5, which is the best-fit point in this scheme, 
XsLMA = 166. 

To summarize, a good fit {x^/doi ~ 1) may be obtained for the atmospheric, LSND 
and CHOOZ data through scheme II. We would like to stress that the best-fit region has 
large LSND angles sin^(26'LSND) — 0.05. Therefore the LSND mass squared difference 
might have observable effects in future measurements, as we would like to discuss next. 



3 Distinguishing between the schemes 

Although comparing the x^'s in the two schemes in order to determine which is the better 
one is not valid, we note that the individual x^'s per degree of freedom are near 1 and 
hence the data seem to be fitted reasonably well. That these x^/dof are similar indicates 
that the goodness of fit for the two schemes are similar, and by itself provides no grounds 
for preferring one scheme over the other. The reason why most authors choose to disregard 
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the LSND results is of course because it still needs confirmation. Since the current data 
do not seem to prefer either of the two schemes, let us examine whether the future data 
will help to shed some light on the identification of the actual mass pattern. 



3.1 Electron neutrino excess in the atmospheric data 

The muon neutrino deficit observed in the atmospheric neutrino experiments can be ex- 
plained well by both the schemes that we are considering. In fact the amount of deficit 
observed is the dominant factor in the determination of the parameters in the fit. Al- 
though we might note that, in the region favoured in scheme II, the LSND angle also 
contributes to a lowering of the down-going muon neutrinos in the multi-GeV range. In 
fact small deficits are indicated by the data and because of the small stastistical errors 
on the multi-GeV muon neutrino ratios the effect contributes to the preference for a large 
LSND angle. Let us now examine the predictions of this scheme about a z/g excess or 
deficit. The data show a small excess, which is accounted for by the SK collaboration by 
an overall flux normalization. The electron excess may be defined as 

ie = ^-l = Pee + rPe^.-l , (16) 

e 



where (N^) is the number of electron neutrinos in the presence (absence) of oscillations, 
and r is the ratio of the original muon and electron neutrino fluxes. The value of r 
depends on the energy as well as on the zenith angle. For the sub-GeV neutrinos, r f» 2 
independent of the zenith angle, whereas for multi-GeV neutrinos, r f» 2 for the near- 
horizontal direction and r «i 3 for the vertical (either up- or down-going) neutrinos. 

In scheme I, there are two factors that determine the extent of excess: Amg and 
S13. The excess due to a non-zero value of Ami, may be written in the form 



Ce(l) = ^Am2(rC23 - 1) , (17) 

whereas the excess due to non-zero S13 may be approximated as ||^| 

ee(2)=P.,3(^4-l)- (18) 

The value of 623 affects ^e(i) and ^e(2) in opposite directions: with an increasing value of 
023, the excess due to the first (second) term decreases (increases). Let us look at the two 
terms separately. 

(i) If 9i3 ^ 0, the solar mass squared difference will produce an excess (deficit) for low 
(high) values of S23. For the sub-GeV data this crossing occurs at s^^ = 0.5. Furthermore 
the excess has the following energy and zenith angle dependence [^: For the sub-GeV 
ratios the excess can be as large as 10-12% for large values of Amg with a positive up/down 
asymmetry but weakly depending on the zenith angle. For the multi-GeV ratios the excess 
is small (5-7 times smaller than for sub-GeV) and the up/down asymmetry is positive. 
The excess is negligible for Am^ < 4 x 10~^. 

(ii) If ArriQ 0, then a non-zero value of S13 will result in an excess (deficit) for high 
(low) values of S23. For multi-GeV up- and down-going neutrinos, this zero crossing point 
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occurs at S23 ~ 0.3. The energy and zenith angle dependence can be described as [0: 
Only in the multi-GeV range and for up-going neutrinos can the excess be substantial 
and hence the up/down asymmetry is positive. 

For both effects we therefore have that for down-going neutrinos in the multi-GeV 
range the electron excess vanishes. Also an excess is accompanied with a positive up/down 
asymmetry. Deficits, though, could give negative up/down asymmetries, but this is dis- 
favoured by the data. Note that for high values of Am|, the best-fit point were found 
to be near S13 ^ and S23 — 0.4, whereas for small Amg the best fit was for non-zero 
(albeit small) values of S13 and S23 — 0.6 (Sec. pTll ). The data therefore seem to prefer 
a small excess in both the sub-GeV and multi-GeV ranges. The ratios for these best-fit 
points are plotted in Fig.^. 

Large values of AttIq would tend to produce a deficit at high S23 for sub-GeV neutrinos. 
Since the data actually show a slight excess, high values of S23 are not favoured with large 
Attt-q. This may be noted from Fig. ^(a), which shows that the SK atmospheric and 
CHOOZ data disallow nearly all the "dark side" (5^3 > 0.5) for Am| = 2 x lO'l On the 
other hand. Fig. ^(b) shows that the observed excess goes on to allow larger values of S13 
for small solar mass differences on the dark side, where the produced excess is positive. 
The present data are however insufficient for us to be able to detect any sign of a possible 
non-zero value of either S13 or Amg. 

In scheme II, the oscillations due to Am^^^^ are averaged out and the electron excess 
becomes 

^ 2sl^cl^{rsl, - 1) + 2rcl,suci2 sin(2^i2 + 2^23) sin^ , (19) 

there are no matter effects and the angles and Am^'s are the same as their vacuum values. 
We have neglected the terms that are more than second order in the small quantity C13. 
For the upper right corner of Fig.^, we have sin(26'i2 + 2^23) — sin(26'atm) — —1 and for 
the lower left corner we have sin(26'i2 + 26'23) ~ 1, so the sign of the second term in (|19|) 
can change. Also, the sign of the first term can change, depending on the value of S23 and 
r, determined by energy range and zenith angle. Small values of S23 are correlated with 
small values of the LSND angle, although the exact value also depends on S13 (533 < 0.4 
is disallowed as seen in Fig. |^ since the LSND angle becomes too small). The best-fit 
points found in Sec. P]2|had S23 ~ 1, and the values with S23 > 0.5 are favoured in the fit. 



Assuming we are in the area with S23 > 0.5, the two terms thus contribute to the 
electron ratio in opposite directions. The first term in (^) is identical for up- and down- 
going neutrinos, when neglecting the small asymmetry caused by the magnetic field of the 
Earth, and produces an excess. The second term vanishes for down-going neutrinos due 
to small L, whereas it gives a small negative contribution for up-going neutrinos. The 
up/down asymmetry is then negative. Furthermore the magnitude of the excess is larger 
for multi-GeV neutrinos, since the value of r is larger. The magnitude depends on the 
interplay between the LSND and the CHOOZ angles. The maximum for a given CHOOZ 
angle is (r — 1) sin^(2^cHOOz)/2. For sin^(26'cHooz) = 0.08 the maximum excess is around 
4% in the sub-GeV range and around 6% in the multi-GeV range. For sin^(26'cHooz) = 
0.03 the effect is only substantial in the multi-GeV range and is around 2%. The sign of 
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the up/down asymmetry is opposite to that obtained in scheme I for both sub-GeV and 
multi-GeV (disregarding the possibihty of a negative asymmetry along with a deficit in 
scheme I ). It should be noted that for small values of the LSND angle a small positive 
up/down asymmetry, in both sub-GeV and multi-GeV ranges, can be obtained in scheme 
II, again normally accompanied by deficits. Nevertheless for some points within the lower 
left corner having S23 — 0.5, a positive up/down asymmetry is obtained along with an 
excess in both sub- and multi-GeV ranges. Hence unfortunately the predictions about 
the electron ratios are not unique is scheme II. 

In the preferred region of scheme II the largest excess is expected to be for down- 
going multi-GeV neutrinos (see Fig. whereas in scheme I this ratio is oneQ Let us 
briefly mention that the large value of r is the reason why the muon neutrino ratios 
are affected less than the electron neutrino ones, in the former the Pufj,ue probability is 
multiplied by 1/r. A clear experimental signal for an LSND mass squared difference is 
therefore an excess of down-going multi-GeV electron neutrinos along with a small muon 
neutrino deficit. With the present statistics, however, it is not possible to distinguish 
between the two schemes. Nevertheless, with the distinct pattern described above, it will 
be made possible by future precise measurements and improvements of the atmospheric 
flux calculation. 



3.2 Long baseline experiments 

Let us look at the possibility of distinguishing between the two schemes in long baseline 
experiments, where the initial neutrino spectrum and flux are better known and hence 
one may expect to have a better handle on the mixing parameters. 

The K2K experiment in Japan started to report their first results . The almost pure 
Ufj, beam (98.2% z/^, 1.3% and 0.5 17^) travels 250 km from the KEK laboratory to the SK 
detector. The average neutrino energy is 1.3 GeV. Hence the value of sin^(Am^^j^L/4£^) 
is large and in scheme II it will dominate over the LSND mass squared terms, which are 
accompanied by small angles. Also, the effect of Am'^^Ya.r negligible in scheme I. This 
already tells us that the predictions for the two schemes will be nearly the same. In 
the following we will neglect the small contamination of electron neutrinos and antimuon 
neutrinos. The expected muon neutrino spectrum at SK, , has been reported in |B2 



after 10^° protons on target (pot). From this spectrum we calculate the number of muon 
and electron neutrino events by performing the integral 

Na = J dE,a^{E,)P^^{E,)^,^{E,) , (20) 

where the probabilities have been calculated using a constant matter density of 2.7 g/cw? 
in the Earth's crust. First we compute the total number of Ve and z/^ events for all points 
within 99% CL to see if it can reveal a difference. The maximum number of Ve events are 



^ The new calculation of the atmospheric neutrinos fluxes indicates that the fluxes are slightly 
smaller than used by SK and the discrepancy is larger for higher energies. This would imply a larger 
electron excess in the multi-GeV range and would therefore seem to prefer scheme II. It could also result, 
though, in an approach to one of the down-going muon neutrinos ratio in the multi-GcV range. 
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found to be 25 for scheme I and 45 for scheme II. However already for sin^(26'LSND) = 0.03 
the maximum number is the same. Therefore in order to see the difference the energy 
spectrum must be measured with a high precision and furthermore the LSND angle must 
be large. Figure |^ shows the expected number of i^e and z/^ events as a function of 
energy for different points within the two schemes, all within 90% CL. The short-dashed 
curve is within scheme II, and chosen to make the number of v^. events small. The long- 
dashed curve, also within scheme II, is chosen so as to make the number of v^. events 
large (equivalently a small number of events). Also two points within the conventional 
case (scheme I) are shown, again one giving a large number of v^. events (dotted curve), 
and one with a small number of v^. events (dot-dashed). The number of events in the 
conventional case falls rapidly above 1 GeV, which is not always the case in the LSND 
scheme. Hence a signal for scheme II could be an excess in the high energy bins. However 
it is clear that the sensitivity in the K2K experiment is not sufficient to either confirm or 
exclude one of the schemes. 

In the planned CERN to Gran Sasso (CNGS) ||31[] long baseline experiment, the main 
purpose will be to detect the appearance of — > Vt. If experimental evidence of a 
appearance is found, a new and very important step in neutrino oscillation experiments 
will be made. A nearly pure beam will travel 732 km from CERN to the ICARUS 
Q and OPERA |Q experiments at Gran Sasso, with a mean energy of 17 GeV. The 
mean energy of neutrinos chosen in this experiment is high so as to cross the r production 
threshold and have a sufficient r production cross section. This large value of the energy 
also turns out to be suitable for distinguishing between the two schemes, since the value of 
sin^(Am^L/4i?) is very small for Am^^j^ or AttIq, but is oscillating for Am'j^g^jj. Hence 
if the LSND angle is fairly large it can result in substantial contributions. 



We again calculate the total number of events in both and channels after 2 • 10 



25 



pot (corresponding to roughly five real years) and per kilo-ton according to (|20|), for all 
points within 99% CL. The measurement of z/^ (see Fig. ^) will not reveal any difference 
between the two schemes unless Aml^^ is well known. Or, turned around, the accuracy 
of Am^^j^ obtained in the case of scheme II is much weaker than that obtained in the case 
of scheme I. 

A large LSND angle will result in more appearance. For OPERA the expected 
sensitivity in the case of a negative search is sin^ 26'^e ~ 1-5 x 10~^ for large Am^. 
The proposed sensitivity of the ICARUS detectors for large Am^ is sin^(2^) > 2.7 x 10^^. 
This is close to the limit obtained in Sec. |2.2| thereby testing nearly the whole LSND region, 
if the decision to build the detector is made. 

The total number of electron neutrinos expected at CNGS as a function of sf^ for 
both schemes is shown in Fig. |^. As can be seen, the conventional scheme I predicts an 
upper bound of < 50 on the number of electron events observed, whereas in scheme II 
the number of electron events can go as high as 400 for sin^(26'LSND) = 0.08. Scheme II 
can thus be identified if the number of events is observed to be large, although it cannot 
be excluded on the basis of a low number of events observed. It is evident that the two 
schemes will have quite different predictions, as the LSND scheme will predict, in most of 
the allowed parameter space, a larger number of electron events than the conventional one. 
Also the energy spectrum of u^. events can be seen in Fig. [TO, using the same four points 
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as in the K2K case above are shown. The energy dependence is also quite different. Hence 
we find that the CNGS experiments should be able to characterize the mass pattern. 

The MINOS experiment |30| will have features quite similar to those of the CNGS 
experiments. The baseline is the same, but the mean energy is lower (around 2 GeV). 
Therefore it will be harder to see any difference between the two schemes, again because of 
the dominance of the atmospheric mass squared oscillation. The total number of electron 
neutrino events is expected to be similar, with a maximum of 20 for scheme I and 28 for 



scheme II (calculated per kilo-ton and per year). In Fig. 11 we plot the energy spectrum 
of the electron neutrino events with a pure i/^ beam. It is seen that only in the case 
of a very large LSND angle will the experiment be able to get signs for an LSND mass 
squared difference. In order for MINOS to be able to distinguish between the schemes, 
the medium or high energy option needs to be employed. Nevertheless the fact that 
the oscillation driven by the atmospheric mass squared difference is dominant allows an 
accurate determination of Am^^^j, contrary to the CNGS experiments. 



4 Summary 

We explore the three-neutrino mixing scheme for solving the atmospheric and LSND 
anomalies, taking into account the constraints from CHOOZ. If the solar neutrino anomaly 
can be accounted for by some exotic mechanism, this scheme can explain all the observed 
neutrino experiments. 

In order to check how well the atmospheric and LSND data can be explained by a 
three-neutrino mixing scheme, we construct a function that takes into account the sub- 
GeV and multi-GeV electron and muon events observed at SK, the 77^ z/g conversion 
probability observed at LSND, and the survival probability observed at CHOOZ. We 
obtain a good fit, with x^/dof ~ 1. In order to compare this goodness of fit with that of 
the conventional fit to the solar, atmospheric and CHOOZ data, we also construct a 
function for the conventional scheme (I). We reproduce the features of the standard three- 
neutrino fits, with x^/dof ~ 1. The two fits are thus equally good: the "goodness of fit" 
criterion, as quantified by x^/dof, by itself does not favour either of these two schemes. 
It is therefore necessary to investigate the often ignored scheme (II) which accounts for 
the atmospheric, LSND and CHOOZ data with three-neutrino oscillations. 

We note some salient features of scheme II. The three-neutrino oscillation does provide 
a modest improvement of the fit with respect to the two-neutrino schemes. Large values of 
the LSND angles are favoured, with sin^(26'LSND) — 0.05, even when including the Bugey 
data. There are almost no matter effects in this scheme, since the values of the Am^'s 
are too large for the Earth's densities to have any effect. 

The two schemes differ in certain respects, which may be exploited for distinguishing 
between them. The ratios observed in the atmospheric neutrinos, for example, should 
display different behaviour. A clear signal for scheme II would be an excess of down-going 
multi-GeV electron neutrinos, accompanied by a small deficit for the down-going multi- 
GeV muon neutrinos. Also the sign of the up/down asymmetry for electron neutrino 
ratios is negative in the favoured region of parameter space, whereas it is positive in the 
conventional scheme. The current data, however, are insufficient to pick out one scheme 
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over the other. 

We also investigate the capabihty of the long baseline experiments — K2K, MINOS, 
CNGS — to distinguish between the two schemes. We compute the Vf. and spectra at 
these three experiments and find that for K2K and MINOS, the LjE value is too large and 
the statistics too small to observe any appreciable difference. However, the final spectra 
at CNGS predicted by these two schemes are very different. We find that the observation 
of the final spectrum, as well as the appearance events that may be observed at 
CNGS, will not give us much useful information regarding the choice of the scheme. The 
apparance of z/g, however, can show some differences in principle, being large in scheme 
II. The observation of a large number of v^. events can rule out the conventional scheme, 
although a low number of events cannot exclude scheme II. 

Scheme II often gets a negatively biased treatment compared with the conventional 
one, mainly because it tries to explain the LSND results, which are not confirmed yet 
by any other experiment. As we have shown in this paper, the goodness of fit cannot be 
reason to prefer one scheme over the other. The ultimate arbiter of the issue is of course 
the experiment. The BooNE experiment will be testing the whole LSND region. With 
MiniBooNE proposed to be fully operative by fall 2002, it will be the first to confirm 
or refute the LSND results. If the LSND result is confirmed, the other experimental 
observations will have to be interpreted in terms of scheme II. 
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Figure 1: The allowed region in the conventional case, (a) for the parameters Aml^^ and 
4:U^3U^3 sin2(2eatm) with 90% (99%) CL marked with a • (•); (b) for the parameters 
sj^ and sig, with 90% (99%) CL marked by a + (*). We have used CHOOZ and SK 
atmospheric data and fitted with four parameters. 
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Figure 2: Allowed region in Aml^^^ and S23 for the conventional case. In both figures the 
+ {^) indicates points within 99% (90%) CL. (a) for all other parameters unconstrained; 
(b) with Am|]^ = 2 x 10~^. We have used CHOOZ and SK atmospheric data and fitted 
with four parameters. 
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Figure 3: The allowed region in sfg and s^^ for scheme I for AttIq constrained to (a) 2 x 
10"^ and (b) 2 x 10~^. We have used CHOOZ and SK atmospheric data and fitted with 
four parameters. 
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Figure 4: Allowed region for scheme II in (a) Aml^^ and — sin^(26'atm) and (b) 

- Uesf sin2(2^cHOOz) and ^U^^U^ ~ sm\2eLSND). Points within 90% (99%) 
CL are marked by a • (■). We have used CHOOZ and SK atmospheric data and fitted 
with four parameters. 
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Figure 5: Allowed region in scheme II for parameters 



and 



'23- 



Points within 68% 



(90%, 99%) CL are marked by • (*,+). We have included the LSND, CHOOZ and SK 
atmospheric data and fitted with four parameters. 
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Figure 6: Ratios as a function of the zenith angle for three best-fit points for (a) sub-GeV 
and (b) multi-GeV. The long-dashed curve is for scheme I constrained to AitIq = 2 x 10~^. 
The solid curve is for scheme I constrained to Amg = 2 x 10~^. The short-dashed curve 
is for scheme II. Upper curves are electron neutrino ratios and lower curves are muon 
neutrino ratios. 
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Figure 7: The expected number of electron and muon events in the K2K experiment as 
a function of the neutrino energy. The solid curve corresponds to no oscillations in the 
muon event plot. The plot is for 10^° pot or roughly 3 years. 
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Figure 8: The total number of z/^ CC events in CNGS calculated for all points within 99% 
CL as a function of Am^^^^ (a) for scheme I, (b) scheme II. The light gray/green points 
in (b) are with sin^(26'LSND) < 0.02 and the dark gray/red points are with sin^(26'LSND) < 
0.08. The plot is for 2 x 10^^ pot or 5 years. 
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Figure 9: The total number of CC events in CNGS calculated for all points within 99% 
CL as a function sfg (a) for scheme I, (b) for scheme II. The light gray/green points in (b) 
are with sin^(2^LSND) < 0.02 and the dark gray/red points are with sin^(26'LSND) < 0.08. 
The plot is for 2 x 10^^ pot or roughly 5 years. 
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Figure 10: The expected number of electron CC events in the CERN to GRAN Sasso 
experiment as a function of the neutrino energy. The plot is for 2 x 10^^ pot or roughly 
5 years. 
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Figure 11: The expected number of electron CC events in the MINOS experiment as a 
function of the neutrino energy. 
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